m 

o 



Oh' 
< 

in 



> 

o 



o 
en 



Vibrational resonance in groundwater-dependent plant 

ecosystems 

C Jeevarathinam'^, S. Rajasekar'^, Miguel A.F. Sanjuan^ 

"'School of Physics, Bharathidasan University, Tiruchirappalli 620 024, Tamilnadu, India 
b 



^\l ' Nonlinear Dynamics, Chaos and Complex Systems Group, Departamento de Fisica, 



Universidad Rey Juan Carlos, Tulipdn s/n, 28933 Mostoles, Madrid, Spain 

Abstract 

We report the phenomenon of vibrational resonance in a single species and a 



Q 
U 

rH ■ two species models of groundwater-dependent plant ecosystems with a bihar- 

monic oscillation (with two widely different frequencies u and Q, Q ^ u) oi 
the water table depth. In these two systems, the response amplitude of the 
species biomass shows multiple resonances with different mechanisms. The 
resonance occurs at both low- and high-frequencies of the biharmonic force. 
In the single species bistable system, the resonance occurs at discrete val- 
ues of the amplitude g of the high-frequency component of the water table. 
Furthermore, the best synchronization of biomass and its carrying capac- 
k> . ity with the biharmonic force occurs at the resonance. In the two species 

c^ i excitable and time-delay model, the response amplitude (Q) profile shows 

several plateau regions of resonance, where the period of evolution of the 
species biomass remains the same and the value of Q is inversely propor- 
tional to it. The response amplitude is highly sensitive to the time-delay 
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parameter r and shows two distinct sequences of resonance intervals with a 

decreasing amphtude with r. 
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1. Introduction 

The study of the response of a system to a small variation of environmen- 
tal changes is important since environmental drivers often fluctuate. The 
fluctuation can be periodic or nonperiodic (noise). In certain ecosystems, 
a small change in one or more environmental parameters leads to consider- 
able changes on their structure and function. It has been noted that many 
systems have relatively high levels of diversity for an intermediate level of a 
disturbance (Roxburgh et al., 2004). The impact of the environmental varia- 
tions/fluctuations has been analysed on food web stability (Vasseur and Fox, 
2007), species coexistence in Holt-McPeek systems (Lai and Liu, 2005) and 
the stability of recovery (Steneck et al., 2002). Without invoking interaction 
between environmental noise and competition, it has been shown that envi- 
ronmental fluctuations enhance coexistence of species which either prefer or 
tolerate different environmental conditions (D'Odorico et al., 2008). 

In the present work, we consider groundwater-dependent plant ecosys- 
tems. It is important to analyse the influence of the variation of various 
environmental factors, particularly, the changes in the water table depth, in 
order to get a deep understanding of the response of various ecosystems. We 



point out that vegetation-water table interactions is very common in many 
ecosystems like wetlands, salt marshes and riparian forests. As a matter of 
fact, it is considered as one of the key mechanisms influencing the dynamics of 
vegetation (Naumburg et al., 2005; Elmore et al., 2006; Munoz-Reinoso and 
de Castro, 2005). Appropriate theoretical models are of great use for explor- 
ing various possible dynamics that can emerge from vegetation-water table 
interactions. In this connection, (Ridolfi et al., 2006, 2007) have proposed 
two vegetation-water table models based on realistic ecological assumptions. 
The first model describes the vegetation biomass dynamics of only one species 
(dominant species). In this model the rate of change of the species biomass 
depends on the existing biomass and the carrying capacity of the system. 
The resultant model is a first-order nonlinear ordinary differential equation 
with a periodic driver. It accounts for multistable states in the dynamics of 
wetland forests and riparian ecosystems (Scheffer et al., 2001). The second 
model describes the two phreatophyte species interacting with a water table. 
In these two models, one of the factors that can change the carrying capacity 
of biomass is the depth of the water table. Change in the water depth due 
to seasonal rainfall oscillations and other sources is represented by a periodic 
function of time. It is also found to display coexistence of two species and 
chaotic dynamics (Ridolfi et al., 2007). 

The infiuence of the environmental variability, treated as a disturbance 
or a kind of noise, in the above two models has been analyzed recently by 
Borgogno and his co-workers (Borgogno et al., 2012). Specifically, they 
have shown the occurrence of stochastic and coherence resonances. When 
a bistable or an excitable system driven by a weak periodic force is subjected 



to an additive noise, it can exhibit an enhanced response at an optimal 
noise intensity. This phenomenon is termed as stochastic resonance (Gam- 
maitoni et al., 1998; McDonnell et al., 2008). Very recently, noise-induced 
spatio-temporal patterns in wetland vegetation dynamics have been reported 
(Scarsoglio et al., 2012). In a subthreshold excitable system, a noise-induced 
resonance can be realized in the absence of external periodic driving and is 
known as coherence resonance (Pikovsky and Kurths, 1997). Interestingly, 
it has been shown that deterministic resonances can be observed even in 
monostable nonlinear systems driven by a biharmonic force in the absence 
of external noise and is called vibrational resonance (Landa and McClin- 
tock, 2000). The analysis of vibrational resonance has received a great deal 
of attention in recent years. Particularly, its occurrence has been investi- 
gated in a spatially extended system in the presence of noise (Zaikin et al., 
2002), Duffing oscillator (Blekhman and Landa, 2004), two-coupled over- 
damped anharmonic oscillators (Gandhimathi et al., 2006) and monostable 
systems (Jeyakumari et al., 2009). Experimental evidence of vibrational res- 
onance was demonstrated in analog simulations of the overdamped Duffing 
oscillator (Baltanas et al., 2003), a bistable optical cavity laser (Chizhevsky 
and Giacomelli, 2006) and an excitable electronic circuit with Chua's diode 
(Ullner et al., 2003). The influence of time-delayed feedback on vibrational 
resonance was studied numerically (Yang and Liu, 2010) and theoretically 
(Jeevarathinam et al., 2011). Further, biharmonic force induced enhanced 
signal propagation was found to occur in one-way coupled systems (Yao and 
Zhan, 2010) and in a coupled network of excitable neuronal systems (Yu et 
al., 2011). 



In the present work, we consider the two groundwater-dependent plant 
ecosystem models of (Ridolfi et al., 2006, 2007) and investigate the emer- 
gence of vibrational resonance. It has been pointed out that the dynamics 
of vegetation have a time scale greater than one season and much greater 
than man-induced periodic disturbances (Ridolfi et al., 2007). We wish to 
mention that high-frequency oscillation of beach water table due to wave 
runup and rundown has been observed and analysed (Waddell, 1976; Li et 
al., 1997). Interestingly, similar high-frequency oscillation of underground 
water table (in addition to the low-frequency periodic oscillation of water 
table due to seasonal variation) can occur due to evaporation, inflow and 
outflow of water and temperature fluctuation. It can also be artiflcially real- 
ized through irrigation or pumping from the aquifer. Furthermore, a water 
table rise and drop can be induced by vegetation removal and planting re- 
spectively. Thus, planting additionally short-lived species interacting weakly 
with species A can also lead to a high-frequency variation of the water depth. 
Therefore, it is realistic to include a biharmonic force in the water table with 
two well-separated frequencies. The flrst model describing the dynamics of 
the biomass V oi a single species has bistable states. When the biharmonic 
force is included in the water depth the system shows an oscillatory variation 
of V. As the amplitude of the high-frequency force is varied, the system ex- 
hibits multiple vibrational resonance with a decreasing response amplitude 
at successive resonances for certain range of flxed values of amplitude of low- 
frequency force. The second model describes the dynamics of two species, 
say A and B, interacting with the water table, and where the evolution of 
B depends on A{t — t) and r is the time-delay parameter. Unlike the single 



species model, the two species model is an excitable system (such a system 
have only one stable equilibrium state, but external perturbations above a 
certain threshold can induce large excursions in phase space, which takes 
the form of spikes or pulses). For a fixed time-delay, both A and B display 
a certain number of resonances when the amplitude of the high-frequency 
force is varied. The resonance profiles of A and B are similar except that at 
resonance the amplitude of A is always much higher than that of B. Here 
the resonance intervals are not sharp but wide. The response amplitude is 
inversely proportional to the period of the variation of A and B. The delay 
parameter r has a strong infiuence on the response amplitude. The response 
amplitude at successive alternative resonances decreases when the value of 
the delay parameter increases. 

2. Vibrational resonance in a single species model 

To start with, first we briefiy introduce the model (Ridolfi et al., 2006, 
2007) in order to prepare the readers for the study of vibrational resonance. 

2.1. Description of the model 

The dynamics of phreatophyte biomass F of a single species (or total 
plant biomass neglecting interspecies interactions) is expressed as (Ridolfi et 
al., 2006, 2007) 

^ = V(Kc-V^), (1) 

where the growth rate of V is assumed to be proportional to the existing 
biomass and the available resources Kc — V^ with V"cc being the carrying 
capacity of the ecosystem, that is, the maximum amount of vegetation sus- 
tainable with the available resources. Based on experimental evidences, an 
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appropriate form of Kc shows a quadratic dependence on the water table 
depth d. Taking into the effect of periodic oscillations in the rainfall regions 
leading to periodic variations of water table depth, (Borgogno et al., 2012) 
considered the form of V^c as 

_ j a [d{t) - dint] [dsnp - d{t)] , if rfjnf <d < d^^p 
I 0, otherwise. 

The form of Kc given by Eq. (|2]) corresponds to the case of phreatophyte 
vegetation that depends on water uptake from the groundwater. In Eq. ([2]) 
d{t) is the water table depth, a is the sensitivity of carrying capacity to 
changes in the water table depth, djnf is the threshold of vegetation tolerance 
to shallow water tables and insufficient aeration of the root zone and dsup 
is the threshold of water depth below which tap-roots cannot extract water. 
The water table depth is given by 

d{t) = do + f3V + F{t), (3) 

where do is the water depth in the absence of vegetation, (3 is the sensitivity 
of the water table to the presence of vegetation and F{t) describes the oscilla- 
tory variation of the water table. The choice F{t) = f cos ut is considered in 
(Borgogno et al., 2012). In the present work we choose F{t) as a biharmonic 
force with two widely differing frequencies: 

F(t) = f cosut + gcosQt, fi > w. (4) 

The explicit time- dependent variation of the water depth can be natural due 
to seasonal rainfall oscillations or man-induced perturbations (pumping from 
an aquifer). 
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Figure 1: (a) The potential U{V) of the system ([T]) for / = 0, g ~ 0. The values of 
the parameters are dg — 0.5m, /3 — 0.5m, a — 26m^^, dmi = 0.6m and dgup = 0.9m the 
values used in (Borgogno et al., 2012). (b) Variation of V"cc with V (thick curve). The 
equilibrium points are the intersections of the bisector (thin straight-line) with the 14c 
curve. The locations of the equilibrium points are marked by solid circles. 



The potential U{V) defined through dV/dt = —dU/dV in the absence 
of F(t) is depicted in Fig. [1^ where do = 0.5m, /3 = 0.5m, a = 26m~^, 
(iinf = 0.6m and dsup = 0.9m the values of the parameters used in (Borgogno 
et al., 2012). U{V) is of a double-well form. The equilibrium states can be 
obtained by setting dV/dt = 0. Vq* = 0, representing the unvegetated state, 
is an equilibrium point. The other equilibrium states correspond to V^c = V. 
In the plot between V^c versus V the intersection points of the line V^c = V 
with the curve of Kc are the equilibrium states. Figure [T]d shows Kc versus 
V. The equilibrium states are Vq* = 0, V* = 0.28526 and V; = 0.56090. V^ 
is the local maximum of U{V) and is an unstable state. Vq and V* are two 
local minima of U{V) and are stable states. When the periodic function F{t) 



is taken into consideration, then the effect of F{t) is to periodically modulate 
the potential U{V). 

2.2. Multiple resonance induced by the biharmonic force 

In the F{t) given by Eq. (|3]), we have assumed that Q ^ u. In this case 
due to the difference in time scales of the low- frequency oscillation / cos ut 
and the high-frequency oscillation gcosQt, the solution of the system ([1]) 
consists of a slow variation of V{t) denoted by V^iow(^) and a fast variation 
Viastity^t). We denote Q^^ and Qq as the response amplitudes of V at the 
frequencies u and Q respectively. A theoretical approach has been developed 
to obtain an analytical expression for Q for certain class of oscillators (Landa 
and McClintock, 2000; Blekhman and Landa, 2004). The threshold existing 
between the carrying capacity V^^ and d in Eq. ([2]) makes it difficult to use 
the theoretical approach to investigate the vibrational resonance. Therefore, 
we compute both Q^ and Qn from the numerical solution of the Eq. ([1]). 
From V(t) the sine and cosine components Q^.s and (5w,c are computed from 
the equations 

(5) 
where T = 2tt/u and n is taken as, say, 200. Then 
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Q^ = ^Qls + Qlc/f- (7) 

Similarly, we can compute Qn from the numerical solution V{t). 

The value of / is significant in observing vibrational resonance in the 
system ([T]). For a fixed value of u and g = 0, small amplitude oscillations 
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of V occur about the coexisting stable equilibrium points V^* and V* for 
I/I ^ 1. Then for fixed values of / with |/| ^ 1 and Q ^ u, when g 
is varied both Q^^ and Qq display one or more resonances depending upon 
the values of the parameters of the system. For g > 0.5, both Q^ and 
Qn ~ and there is no further resonance. Therefore, we present the results 
ioY < g < 0.5. Figure [2] shows the dependence of Qi_j and Q^ for three 
fixed values of / with a; = 1 and i7 = 10. In all the cases a nonmonotonic 
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Figure 2: Response amplitudes (a) Qui and (b) Qn of plant biomass V for three fixed 
values of the parameter / with cj = 1 and 17 = 10. The values of the other parameters in 
the system ([1]) are as in Fig. [TJ 
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variation of both Q^ and Q^ occurs. In Fig. [JlQw and Q^ are ^ Qioi g > 0.5. 
For / = 0.08 as g increases from a small value, Q^j initially decreases then 
increases sharply, reaches a maximum at f^ = (7^^^ = 0.069 and then sharply 
decreases to a lower value. As g increases further, Q^^ becomes maximum at 
two other values of g, namely g = g^^^ = 0.126 and g = g^-^^ = 0.312. There 
are three resonances. We note that a resonance at the low-frequency of the 
water table oscillation is induced by an appropriate value of the amplitude 
g of the high-frequency oscillation of the water table. For this reason the 
above resonance phenomenon, that is the occurrence of a maximum of Quj, 
is termed as vibrational resonance. For / = 0.05 also three resonances occur. 
For both / = 0.05 and 0.08 the first resonance is the dominant and the value 
of Q at the successive resonances decreases. On the other hand, the width of 
the successive bell-shaped resonance curves becomes wider and wider. When 
the value of / increases the resonance peaks move towards lower values of g. 
For / values above a critical value, the first two resonance peaks disappear. 
Moreover, the value of Q at gf = is substantially enhanced with an increase 
in the values of / producing a resonance without tuning. 

In the nonlinear oscillators driven additively by a biharmonic force, Qq 
monotonically increases when g increases and there is no resonance-like vari- 
ation of it. We wish to remark that in the system ([1]), the biharmonic force is 
not an additive force but it is in the expression for V^c- Therefore, one wish 
to know the response of the system ([T]) and its relative strength at the high- 
frequency Q of the water table oscillation. For this purpose, we computed 
and presented the variation of Q^ with g. Interestingly, Qq also displays 
resonance. This is shown in Fig. |2]3. The dependence of Q^ on g is similar 
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to that of Qi^. For / = 0.05 and 0.08 there are two resonances with Q^ and 
both occur at the same values of g at which Qw becomes maximum. Qq also 
shows resonance without tuning. Though the resonances associated with the 
two frequencies u and Vt occur at the same value of the control parameter 
g, at resonance Q^ is much higher than Qq. We note that Q^ and Qq, are 
proportional to the Fourier coefficients of the periodic terms with the fre- 
quencies u and VL respectively in the Fourier series of V{t). Since the Fourier 
coefficients decay with increase in the frequency and because Vt^ u the re- 
sponse amplitude Qq is much less than Qt^ for each value of g in Fig. [2l It is 
noteworthy to mention that aX g = g^-^ an enhanced vegetation can be real- 
ized over two time intervals, one with the low-frequency u and another with 
the high-frequency Q. These two frequencies need not be commensurable. 

2. 3. Dynamics of V and V^c at and far from resonances 

The occurrence of resonance in the system ([T]) can be understood by 
analyzing the influence of the control parameter g on the evolution of the 
plant biomass V and the carrying capacity Kc of the ecosystem. Figure [3] 
shows V{t) and the carrying capacity V^c(^) of vegetation biomass for several 
fixed values of g. In the absence of F{t) there are two stable dynamical 
states V^* and V*. The values of / and u are chosen in such a way that in 
the absence of a high-frequency oscillation of F{t) the variable V{t) oscillates 
either about V^* or V* depending upon the initial condition and there is no 
transition between these two states. When g is varied from then for very 
small values two oscillatory states coexist. The equilibrium states about 
which oscillation occurs are perturbed by F{t). For the entire range of values 
of g considered in Fig. [2] the evolution of V is periodic. V{t) is said to be 
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Figure 3: The time evolution of (a)-(f) V{t) and (g)-(l) Vcc{t) of the system ([T]) at few 
values of g with / = 0.08, w = 1 and ft — 10. The values of the other parameters are as 
in Fig. [TJ The dashed curves represent f cos cut. For comparative purpose the amplitude 
of f cosujt is suitably multiplied by a factor. 
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periodic with period t' if V{t + t') = V{t) for some finite and nonzero value 
of t' after leaving initial transient evolution of it, say for example leaving 
V{t) for < t < 100 X T(= 27r/a;). The period of V{t) is found to be 
T = 1t\ juj, the period of F(t) given by Eq. (jl]). In Fig. [3^ for gf = 0.02 far 
before the first resonance, for the initial condition chosen near V^*, though 
both Vif) and /cos cut are periodic with the same period T the forms of 
both of them are different. Over one drive cycle both Vif) and V^c(i^) (see 
Fig. [3^) have more than one dominant maximum while F{f) has only one. 
There is no synchronization between V (as well as Kc) and fcosut. For 
g = g^-^-^ = 0.069 in Figs. [3b and h the center of oscillation is shifted. V and 
Kc are synchronized with fcosut (except that there is a phase difference). 
This feature of V leads to the first resonance. The numerical analysis shows 
an absence of resonance at g = 0.069 for the vegetation dynamics about 
the other low value stable state (Vq). The value of Q of the corresponding 
dynamics is very small. 

Next, at (7 = 0.07 the high-frequency oscillation of the water table induces 
a transition of V about V^* to Vq. When the vegetation biomass is close to 
V^* = 0, it is easy to note from Eqs. ([2]) and ([3]) that, Kc = for most of 
the time over a drive cycle of F{t). Consequently, the dynamics is confined 
near Vq. There is no transition between the two states V^* and V*. This is 
shown in Figs. [St and i ioi g = 0.075. Though both V and V^c are quite 
synchronized with fcosut, there is no resonance because V is trapped into 
the unvegetated state. Because the center of oscillation of V is shifted from 
the vegetated state to the unvegetated state both Q^^ and Qc make a sudden 
jump from a higher value to a smaller value at g = 0.07 (see Fig. [2]). When 
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g increases further, the time intervals in which V^c = decreases. In Fig. [3] 
corresponding to g = 0.126 the total time over which V^c = is ~ T/2 and 
the line joining the maxima of \4c (indicated by solid circles connected by a 
line) varies sinusoidally and one can clearly notice synchronization between 
V^c and /cos cut. A second resonance occurs at this value of g. In Fig. [31i 
V oscillates about the perturbed equilibrium state V*. For g = 0.23 rapid 
oscillations of V and V"cc (Figs. [3^ and k) take place and Q^ ~ 0. At g = 0.312 
a third resonance with Qcj much less than that of the first two resonances 
occurs. However, we can clearly see the synchronization between V and 
f cosut in (Fig. [3]F). 

In this section so far we reported the results for a; = 1, fi = 10 and 
/ = 0.08. We numerically studied the occurrence of resonance for a wide 
range of fixed values oi u, Q and / and thereby varying the control parameter 
g. Figures UK and b show the variation of Q^^ with (7 for w G [0.5, 1.5] where 
i7 = 10 and fi G [5, 15] where a; = 1 respectively. Three resonances occur for 
a wide range of values of u and Q. Figures Ht and d illustrate the influence of 
the parameter / on resonance for two sets of values of u and Q. In both the 
cases we observe one or more resonance for / ^ 1 . The values of g at which 
resonance occur and the corresponding value of Q vary with the parameters 
of the system. 

3. Vibrational resonance in the two species model 

In the previous section, we have described the vibrational resonance asso- 
ciated with the vegetation dynamics of a single species interacting with the 
water table. The present section is devoted to a two phreatophyte species A 
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Figure 4: Variation of Qi^ as a function of (a) a; and g for fi = 10, / = 0.08, (b) ft and 
5 for w = 1 and / = 0.08, (c)-(d) / and g for two sets of fixed values of w and il. The 
values of the other parameters are as in Fig. [TJ 

and B interacting with the water table. 

3.1. Description of the model 

In the two species (denoted by A and B) model proposed by (Ridolfi et 
al., 2006, 2007) the species A is assumed to be dominant over the species B. 
That is, in the absence of interactions with the water table, A tends to its 
maximum density while B tends to disappear. The values of A{t) and B{t) 
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are normalized with respect to their maximum value. Frequently, the logistic 
law is chosen for the growth of A and B. The model system is thus 

— = aAA{V,A-A), (8) 

^ = aBB{V,B-A-B), (9) 

where a^ and a^ are the coefficients determining the response rate of A and 
B respectively and Vca and Vcb are the carrying capacities of the species A 
and B respectively. Vca and Vcb depend on the depth d(t) of the water table. 
Introducing the change of variables A = as/o^A and t = a At' (which makes 
time a dimensionless quantity), dropping the prime in t' and assuming that 
the dynamics of B depends on A{t — T) where r is a time-delay, the Eqs. dHHH} 
become 

A A 

— = A{V^A-A), (10) 

^ = \B{V,B-Ait-T)-B). (11) 



In Eqs. Umlii 

Vci = 9[d — dmm,i] ■ d [(iinax.j — d] , (12) 

where i = A,B, 6 [s] is the Heaviside function, that is, 6* [s] = 1 for s > 
and for s < 0, dmin,i and (imax,j are the minimum and maximum water table 
depths tolerated respectively by the species i with dram,B < dmm,A < '^max,s 
(species A needs a deeper aquifer than species B) and 

d{t) = do + /3aA + /3bB + fcosujt + g cos nt, (13) 

where do is the water table depth in the absence of vegetation {A = B = 
0) and Pa and Pb are coefficients that weight the control exerted by the 
vegetation on the water table depth d. 
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In the above two species model one of the species is assumed to be sub- 
dominant. Further, the feedback of A on the evolution of B is taken as a 
time-delayed term. For a variety of phreatophytes these assumptions can be 
realized. Dominant species A can be regarded as plants with tap-roots able 
to penetrate through relatively deeper than those of the species B. In the 
literature of plants, approximate maximum lengths of tap-roots of phreato- 
phytes are reported. Some of the plants with deep penetrating tap-roots with 
length more than 15m are mequite, camelthorns, grease wood, and purple 
medic. Plants such as black grease wood and banksia (173 species) have 
tap-roots with length in the range of 5- 10m. Examples of plants with tap- 
roots of short length about l-5m are saguaro, creosote bush, ocotillo brittle 
bush, sagebrush, alder and chamisa. One can identify appropriate species 
of types A and B. In this connection we wish to cite that the experimental 
analysis carried out on holm oaks and cork oaks indicated that the higher 
water status leading to more effective drought avoidance of former is due 
to their deeper root systems compared to the latter (David et al., 2007). 
A field experiment was performed on the two phreatophytic plant species 
Alhagi sparsifolia (camelthorns) and Karelinia caspia occurring around the 
river Oasis at the southern fringe of the Taklamakan desert (Vonlanthen et 
al., 2010). Both the species occur at sites with distances to the ground water 
table upto 12m while only Alhagi sparsifolia occurs at distances upto 17m. 

The motivation for introducing the time-delay term A(t — T) in the system 
( 1T01ITT|) is to take into account the fact that the changes in the population 
of a species, generally, will not have immediate effect on the growth of own 
population and on the interacting species. The effect will be realized after 



a time-lag. The effect of time-delay has been studied in population models 
and vegetation dynamics (Kuang, 1993; Wang et al., 2011). Reports on the 
analysis of influence of various factors on the growth of vegetation dynam- 
ics based on AVHRR (Advanced Very High Resolution Radiometer) images 
indicate that the time-lag can be few days to few months (Richard and Poc- 
card, 1998; Li et al., 2002; Wang et al., 2006; Farajzadeh et al, 2011). A 
reasonable value of time-delay can be of the order of l/ag where as is the 
coefficient determining the response rate of the species B. 

3.2. Multiple resonance 

For our numerical study we fix d^ = Im, (imin^A = 1.5m, (imax.A = 2.5m, 
dmin,B = 0.5m, (imax.B = 2m, /J^ = 0.51m, f3B = 0.9m, aA = as = ld"S 
oj = 0.5 and Q = 5. When / = and g = the system (llOtllip has 
two equilibrium states. {A*,B*) = (1,0) is stable while {A*,B*) = (0,1) is 
unstable. {A*, B*) = (1,0) remains as a stable equilibrium point for / < 0.01 
when u = 0.5, g = and r = 1. The system exhibits excitable dynamics for 
/ > 0.01. 

We choose the value of / less than 0.01 so that in the absence of high- 
frequency oscillation of water table the system is in the stable equilibrium 
state. It is important to investigate the response of the system for a wide 
range of values of the control parameters /, co, g, fl and r. Because such 
a study in the five parameters space is time consuming we restrict to {g, r) 
parameters space. We studied the response of the system for a wide range 
of values of g for several fixed values of the parameters of the system. For 
g > 0.4 there are no resonances and the response amplitudes of both the 
species A and B are ~ 0. Therefore, we consider the range of g as [0,0.4]. 
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Before presenting the results in {g, r) parameter space first we discuss the 
presence of resonance for fixed values of /, u, fl and r thereby varying the 
parameter g. 

Figure [5] presents the numerically computed response amplitudes Qa{^) 
and Qb{^) as a function of g for three fixed values of / with r = 1. There 
are few interesting results: 
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Figure 5: Response amplitudes Qa and Qb of species A and B versus the parameter g for 
three fixed values of /. The values of the other parameters are do = Im, dniin,A = 1.5m, 
dmax.yi = 2.5m, dmin.B = 0.5m, dm&yi,B = 2m, Pa = 0.51m, Pb = 0.9m, a^ = as = ld"\ 
ijj = 0.5. il = 5 and r = 1. 
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• The response amplitudes of the species A and B display similar varia- 
tion, however, Qa > Qb, that is, the species B is subdominant. 

• An interesting result is that even though species A is dominant, when 
A exhibits resonance the species B does not disappear but also displays 
a resonance and, moreover, in Fig. [5] we notice that Qb ~ Qa/'2 in the 
resonance region. This is due to the vegetation-water table interaction. 

• Multiple resonance occurs. Qa and Qb are not maximum at discrete 
values of g, but they are almost constant over a range of values of g. 

• Qa (and Qb) values at successive resonances are not equal. In the 
single species model also response amplitude values are not same at 
successive resonances. We wish to point out that in bistable systems 
driven additively by a biharmonic force the response amphtude at suc- 
cessive resonances are found to be the same (Landa and McClintock, 
2000; Rajasekar et al., 2010). 

• Qa (and Qb) values at resonances decrease by increasing the value of 
/, while the width of the resonance interval increases. 

To understand the occurrence of a multiple resonance and almost plateau 
regions of resonance profile, we consider the nature of the evolution of the 
system. The system fllOtilip is a system of two-coupled first-order nonlinear 
differential equations driven by a periodic force. Such a system is capable 
of exhibiting different types of nonlinear dynamics including chaotic dynam- 
ics (a nonperiodic and bounded evolution of a system with high sensitive 
dependence on initial conditions). In the system f fTOlfTTj) in the absence of 
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time-delay and high-frequency oscillation of water table chaotic dynamics is 
found to occur when the amplitude of the low-frequency oscillation of the 
water table is varied (Ridolfi et al., 2007). However, in the system f lTOlfTTl) 
for the parametric choices considered in the present work any route to chaotic 
dynamics is not observed. 

In Fig. [5] the variation of Qa and Qb is shown for g G [0, 0.4] only. For 
g > 0.4, as mentioned earlier, we found Qa and Qb ~ and hence we 
restrict ourselves to the interval < g < 0.4. In this interval of g either 
stable equilibrium state or periodic variation of A and B is found depending 
upon the value oi g. To identify the period of A{t) we collect the values 
of A{t) at t = nT where n = 1, 2, ■ ■ ■ , m{= 500), T = 27i/u and designate 
them as A„. If A{t) is periodic with period T then Ai = A2 = ■ ■ ■ = Am- 
For a period-2T variation of A{t) we observe Ai = A3 = ■ ■ ■ = Am-i = o-i, 
A2 = A4 = ■ ■ ■ = Am = CL2 and ai 7^ 02. Similarly, we can define higher 
periods of A and identify the periodicity of A. In Fig. we have plotted the 
period of A together with Qa for / = 0.005. The period of A is a complicated 
function of g. 

Interestingly, the period remains the same in the regions of g where Q is 
almost a constant. Different plateau regions correspond to different constant 
periods of A. For g < 0.007 we observe {A,B) — )• (1,0) as t increases. The 
species biomass A attains its maximum while that of the species B failed 
to survive. Even though the water table depth d{t) varies periodically, the 
evolution of A and B is not oscillatory but approaches the equilibrium state 
for g < 0.007 as shown in Figs. [7^ and e for g = 0.005. This is because the 
oscillatory variation of d{t) does not strictly appear as an additive driving 
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Figure 6: Plots of (a) the response amplitude Qa and (b) the period of A (in units of T) 
versus g for / = 0.005. In the subplot (a) three most dominant intervals of g with large 
values of Qa are marked by the numbers 1, 2 and 3. The values of the other parameters 
in the system (jlOlllip are as in Fig. [5] 

force in fllOtllip but its influence is on Vca and Vcb- Equations fllOHllI) can be 
regarded as a parametrically driven system. Moreover, Vca and Vcb are not 
sinusoidally varying function of time. They take the values or 1 depending 
upon the values of the various parameters of the system and the value of t. 
The most dominant resonance interval is g & [0.157, 0.21]. In this interval 
of g the period of A is T (see Fig. Eb). Figures [7b and g show the evolution of 
A and B with time for g = 0.2. We can clearly notice the best synchronization 
of A and B with the external drive fcosut. The next dominant resonance 
interval is g E [0.059,0.071]. Figures [6b indicates that the period of A in 
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Figure 7: Temporal behaviour of (a)-(d) A and (e)-(h) B for four fixed values of g with 
/ = 0.005. The values of the other parameters in the system (|10mip are as in Fig. [S] The 
dashed curve represents / cos ujt and its amplitude is suitably magnified and its center of 
oscillation is shifted for illustrative purpose. 

this interval of g is 2T. In Figs. Uh and [Tf for g = 0.065 we observe pulse- 
like solution with period-2T. Period-3T oscillation of A and B occurs in the 
interval g G [0.037, 0.047] (not shown in Fig. [7]) and this region constitutes 
the third dominant resonance region. 

For g > 0.24 the evolution of A and B exhibits a rapid oscillation (as 
shown in Figs. [Tli and h ioT g = 0.242) and both Qa and Qb decays to 
zero implying degrading of the response of the biomasses A and B with 
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increase in the value of g. Outside the three dominant resonance regions, 
we observe few small intervals of g where Q becomes almost constant but 
relatively very small due to higher periodicity of variation of A and B. In 
Fig. [7] we notice that as g increases from a small value the evolution of A 
and B undergoes a transition from equilibrium state — )■ pulse-like solution — >■ 
rapid oscillatory solution. Further, when A becomes maximum (minimum) B 
becomes minimum (maximum). From the above, we point out the following 
features: 

• Q is very small if the species biomass remains a constant (see Figs. [7^ 
and e) or periodic with large period or it oscillates rapidly (see Figs. Eli 
and h). 

• Q is larger, for example the regions 1, 2 and 3 in Fig. |6^, when the 
period of the pulse-like evolution of the biomass species is lower. 

• Q remains almost constant over an interval of g (regions 1, 2 and 3 of 
Fig. [6^) if the period of A remains the same. 

It is noteworthy to compare the mechanisms of vibrational resonance 
in the FitzHugh-Nagumo (FN) equation, a well studied excitable system, 
and the system (IIOHIII) . In the FN equation when the amplitude of the 
high-frequency force varies, a resonance takes place when the waiting time 
(T^), the time the system spends around the equilibrium point between two 
consecutive firing, is T/2 (Ullner et al., 2003; Hu et al., 2012). In the noise- 
induced stochastic resonance also a resonance occurs when T„ ^ T/2. In 
the system (llOtilll) . Q is inversely proportional to the period of evolution of 
biomasses. 
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3.3. Effect of a time-delay on the resonance 

So far, we have focused our analysis for a specific value of the time-delay 
r(= 1). Now, we present the effect of time-delay on the resonance. 

First, in the parameter space {g, r), we identify the regions where Q(t, g) > 
(5(0, g) for both species. We divide the {g, r) parameter space in the interval 
g e [0, 0.4] and r G [0, 12] into 100 x 100 grid points. We collect the grid 
points for which Q{T,g) > Q{0,g). The result is the Fig. [SI Both Figs. [H^ 
and b corresponding to the species A and B respectively are almost similar 
except near g = (far before resonance) and near g = 0.4 (far after the 
last resonance) and in these two regions Qa and Qb are very small. The 
time-delay r has a strong influence on the response amplitudes of A and B. 

Next, Fig. [9] features the colour-contour plots of the dependence of Qa 
and Qb on g and r. Multiple resonance is found to occur for a wide range of 
fixed values of r. Moreover, Fig. [H] clearly demonstrates the strong influence 
of T on the variation of Q. For fixed values of g the response amplitudes Qa 
and Qb are not periodic with r. Figure [10] depicts the dependence of Qa 
on the time-delay r for (7 = 0.1. Qa (as well as Qb) exhibits a sequence of 
resonance intervals. The period of A (as well as B) in the first resonance 
interval is 2T. The periods of A in the other consecutive intervals are 3T, 
4T, ■ ■ ■ respectively. That is, the system exhibits a resonance sequence with 
period adding dynamics. The resonance intervals with period of A being even 
(odd) integer multiples of T are marked by open (solid) circles in Fig. [TOl The 
values of Qa of these two sequence of resonance intervals decrease rapidly 
with r. Further, we observe that Qa{'^) does not varies periodically with 
r, that is (5^(t + «) 7^ Qa(t) for some finite nonzero value of a. However, 
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Figure 8: Regions (marked by black colour) where (a) QA{T,g) > Qa{t = 0,g) and (b) 
QB{T,g) > Qb{t — 0,g) for the system pOlllip . The values of the parameters are as in 
Fig. El 

the resonance intervals marked by solid circles and open circles occur at a 
regular interval of delay-time ~ T. Similar dependence of Q on r is found 
to occur for other fixed values of g. We note that in nonlinear oscillators 
with time-delayed feedback and driven additively by a biharmonic force Q is 
shown to vary periodically (Yang and Liu, 2010; Jeevarathinam et al., 2011). 
In this section we presented our analysis on the system (llOtllip for u = 0.5 
and Q = 10. Results similar to these fixed values of u and Q are observed 
for a wide range of values of u and Q. The number of resonances, the values 
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Figure 9: (Colour in online) Colour-coded dependence o^Qa and Qb of the system (|10j|lip 
on the parameters t and g for / — 0.005. The values of other parameters are as in Fig. [5] 

of Q at resonances and the values of g at which resonances occur depend on 
the values of the other parameters of the system. 

4. Conclusion 

Several studies have reported the occurrence of vibrational resonance in- 
duced by a two-frequency periodic force in physical and biological nonlinear 
systems. The present work is devoted to the analysis of the effect of a bi- 
harmonic type variation of the water table depth in two phreatophyte plant 
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Figure 10: Qa of the system pOJIlip as a function of time-delay r for g = 0.1. The values 
of other parameters are as in Fig. [5l The period of A in units of T(= 27r/w) is marked 
for first few resonance intervals. The resonance intervals with period of A equal to even 
(odd) integer multiples of T are marked by open (solid) circles. 

ecosystems. We have considered a single species and a two species model sys- 
tems. The former has bistability, while the later is an excitable system. Our 
study shows how a very simple deterministic periodic variation of the water 
depth is able to give rise to a great increase in the response of vegetation 
dynamics of single species and two species ecosystems. 

In both overdamped and underdamped nonlinear oscillators exhibiting 
multiple vibrational resonance (Landa and McClintock, 2000; Baltanas et 
al., 2003; Jeyakumari et al., 2009; Yang and Liu, 2010; Jeevarathinam et al., 
2011), at all resonances the long time motion of the systems is found to be 
periodic with period T of the low-frequency force when the ratio Q/u is an 
integer. In the single species system ([T]) the period of evolution is T at all 
the resonances, however, the value of Q at successive resonances decreases. 
In the two species model, the period of the system is different at resonances. 
Further, the response amplitude Q is found to be inversely proportional to 



29 



the period of the evolution of the species biomass. Analysis of vibrational 
resonance in systems described by different kinds of evolution equations can 
lead to a deeper understanding of the phenomenon. 

We believe that analysis of nonlinear phenomena such as chaos, stochas- 
tic and vibrational resonances, in the theoretical models of vegetation-water 
table interactions may motivate the experimentalists to perform experiments 
with controlled variations of the water table depth in a small scale. Such 
studies not only would explore the ways of enhancing the biomass response, 
but also help to improve the theoretical models based on experimental ob- 
servations. 
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